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Abstract 

We show that heavy meson hadronic decay widths satisfy quark-hadron duality when smeared 
over the heavy quark mass, M, to an accuracy of order 1/M^. 
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Quark-hadron duality has been part of the lore of strong interactions for three decades. 
Bloom and Gilman[^, ^ (BG) discovered duality in electron-proton inelastic scattering. 
There, the cross section is given in terms of two Lorentz invariant form factors Wi and 
W 2 which are functions of the invariant mass of the virtual photon, q^, and the energy 
transfer to the electron, ly. Gonsidering the form factors as functions of the scaling variable 
u = q‘^/2Mi>, they compared the scaling regime of large (and large z/) with the region of 
hxed, low q^. They determined that, for each form factor, the low q^ curves oscillate about 
the scaling curve, that identihable nucleon resonances are responsible for these oscillations 
and that the amplitude of a resonant oscillation relative to the scaling curve is independent 
of q^. Moreover, they introduced sum rules whereby integrals of the form factors at low and 
large q^ agree and noticed that the agreement was quite good even when the integration 
involved only a region that spans a few resonances. 

Poggio, Quinn and Weinberg[|] (PQW) applied these ideas to electron-positron annihi¬ 
lation. While BG compared experimental curves among themselves, PQW compared the 
experimental cross section to a scaling curve calculated in QGD. They noticed that the 
weighted average of the cross section ct(s). 


{a{s)) = — ds — - .. ^ (1) 

TT Jo [S — S)^ -|- 

is given in terms of the vacuum polarization of the electromagnetic current with complex 


argument. 


and argued that one can safely use perturbation theory to compute this provided A is large 
enough. This procedure was better understood with the advent of Wilson’s]^ Operator 
Product Expansion (OPE). It is interesting to point out that the prediction of PQW based 
on the two generations of quarks and leptons known at the time did not successfully match 
the experimental results. When PQW allowed for additional matter they found a best 
match if they supplemented the model with a heavy lepton and a charge 1/3 heavy quark, 
anticipating the discoveries of the tan-lepton and b-quark. 

In an attempt to understand the origin of quark-hadron duality we have computed both 
the actual rate and its “scaling limit” from hrst principles in special situations. In Ref. [Q we 
computed the semi-leptonic decay rate and spectrum for a heavy hadron in the small velocity 
(SV) limit. We showed that two channels, B —>■ Dev and B —> D*ev, give the decay rate 


2 



to first two orders in an expansion in l/irtb and that to that order the result is identical to 
the inclusive rate obtained using a heavy quark OPE as introduced in Ref. P . The equality 
holds for the double differential decay rate if it is averaged over a large enough interval of 
hadronic energies. The computation demonstrates explicitly quark-hadron duality in semi- 
leptonic R-meson decays in the SV limit, but really sheds no light into the mechanism for 
duality. In particular, it is puzzling that duality holds even if the rate is dominated by only 
two channels. 

More recently we attempted to verify duality in hadronic heavy meson decays. In Ref. 0 
we considered the width of a heavy meson in a soluble model that in many ways mimics 
the dynamics of QCD, namely an SU{Nc) gauge theory in 1 -|- 1 dimensions in the large 
Nc limit. This model, hrst studied by’t Hooftp], exhibits a rich spectrum with an inhnite 
tower of narrow resonances for each internal quantum number, making the study of duality 
viable. We considered a ‘R-meson’ with a heavy quark Q and a light (anti-)quark q of 
masses Mq and m, respectively, which decays via a weak interaction into light qq mesons. 
To leading order in 1/W the decay rate is dominated by two body hnal states: if 7ij denote 
the tower of gg-mesons, the total width is given by r(R) = I]r(R —>■ TTjTrk), where the 
sum extends over all pairing of mesons such that the sum of their masses does not exceed 
the R mass, + ^k < Mb- The main result of that investigation was that there is rough 
agreement between r(R) and the decay rate of a free heavy quark, r(Q). When considered as 
functions of Mq the quark rate is smooth but the meson rate exhibits sharp peaks whenever 
a threshold for production of a light pair opens up. This is due to the peculiar behavior 
of phase space in 1 -|- 1 dimensions, which is inversely proportional to the momentum of 
the hnal state mesons. Nevertheless, in between such peaks it was found that the relation 
r(R) = r(Q)(l -|- 0.14/Mq), in units of g^Nc/n = 1, holds fairly accurately. 

Recently]^] we considered the effect of local averaging on the results of Ref. 0. The main 
result is that when averaged locally over the heavy mass Mq the agreement between r(R) 
and r(Q) is parametrically improved. In fact, for the averaged widths we found 


(r(R)) ^ (r(Q)) 


0.4 5.5 


( 3 ) 


Remarkably, the correction of order 1/Mq has disappeared. 

In this paper we demonstrate that when averaging over Mq the corrections of order 1/Mq 
are absent. The argument we present is very general and applies both to the’t Hooft model. 
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FIG. 1: Analytic strcuture of the Green function T{Q) = T{Mv, Qv) in the complex Q plane. Guts 
on the real line are depicted by a heavy solid line. 

explaining the nnmerical observations of @ , and the phenomenological relevant case of fonr 
dimensional QCD. The central idea is simple. In a heavy qnark effective theory the fonr 
qnark operator describing a weak i?-meson hadronic decay, is 

where is the heavy qnark held with velocity v. The exponential factor, which acconnts 
for the large momentnm carried by the heavy quark, plays the same role as an insertion of 
external momentum exp(—ig ■ x) with the specihc choice q = Mv. Thus one can use dis¬ 
persion relations to relate the decay amplitude to Green functions with complex momentum 
where an OPE is valid, much like the procedure for semileptonic decays in Ref. p|. The 
resulting relation has then the form of a mass averaged amplitude in terms of a systematic 
OPE. 

Consider the Green function 

T(g,p) = i J d^x e^^-^{B{p)\T{:H\x)n{Q))\B{p)) (5) 

where the B momentum is p = Mv and is the term in the weak Hamiltonian density 
responsible for hadronic B decay: 

AGf 

'X = -^VcbKdClYbdLlF'^L- ( 6 ) 

A simple calculation gives 

ImT = ^7r(27r)35^(g + p-px)KX|R:(0)|H)P (7) 

+ ^7r(27r)^(5^(g -p + Px)K^|7T^(0)|.B)P (8) 

X 
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Hence, the analytic structure of T{Q) = T{Mv, Qv) is as shown in Fig. |l|. The two real axis 
cuts are associated with the two time orderings of Ti and 1-0. The discontinuity across the 
hrst cut, which runs from Q = Qi = to inhnity, is related to the inclusive 

decay rate of the B meson. For Q > ms there is also a contribution from states with two 
B mesons. The discontinuity across the second cut, running from Q = Q 2 = — — m£, to 
negative inhnity, is related to a process with two units of H-number in the hnal state. In 
addition, a pole at —rriB^ is not shown. The decay rate is obtained as the discontinuity at 
g = 0, 

F(H) = —ImT|,=o. (9) 

tub 

While T may be computed perturbatively when the complex momentum is sufficiently 
away from the real cut, the computation of F(H) requires T at one point on the cut itself. 
This has been the main impediment to computing the decay width. In processes such as 
e’''e“ annihilation into hadrons or in semileptonic B decays, an integration over q allows one 
to use a dispersion relation that relates an integral of the discontinuity of T on the real axis 
to the value of T in the complex plane. But in this process q is hxed. 

Our solution to this problem makes use of the observation above that when computing 
T in an effective theory for the static heavy quark the momentum of the heavy quark, Mv, 
and the external insertion of momentum, q, enter all expressions in the precise combination 
q + Mv. Therefore, one may still use a dispersion relation integrating over q, and this will 
have the same effect as an integral over M. The result is a perturbative expression for an 
integral over the mass of the decay width. 

We dehne a Green function in the effective theory similarly, 

f{q,v) = i J d^x 0^^+^^^-^{H{v)\TOfO{x)n{t]))\H{v)) ( 10 ) 

Here H{v) is the state corresponding to the static quark with four velocity n, with a non¬ 
standard normalization (independent of M) as is appropriate in the effective theory. To this 
order, the weak Hamiltonian in the effective theory, 7d, is the weak Hamiltonian B. with the 
quark b replaced by the effective theory static held h^. It follows that 

Imf = Y,'^i‘^T^fS\(l + Mv-px)\{X\n{0)\H)\0 (11) 

A second term, of the form 

Y,M2^fS-'(q + Mv+px)\{X\HimH}f. ( 12 ) 

X 
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FIG. 2: Contour for the integral in the complex Q plane leading to Eq. (13). The thick line on the 
real axis denotes the cut of the Green function T{Q). 


is absent because the intermediate state X is required to have two units of 5-number, which 
is excluded from the effective theory. This also means that the left cut in Fig. is absent 
for T. The effective theory, while unable to properly reproduce the full Green function 
T{p,q), does provide a systematic approximation to physical quantities like the hadronic 
decay width, thus r(5) ImT|q=o- 

To leading order in 1/M one has an additional result: since the mass enters only in the 
combination q + Mv, the decay rate for a heavy meson of heavy quark mass M + 6M is 
r(M -I- 5M) lmT\q=sMv If we dehne f{Q) = T{Qv,v), which depends implicitly on M, 
then to leading order in M the dependence on M and Q is only through the combination 
M + Q. It is straightforward now to use standard methods of analysis to relate an integral 
of T to the Green function of complex argument. Using the contour in Fig. |], we have 


d 


'n—l 


f(z) 


{n — 1)! dz^~^ {z + iXy 


+ 


d 


'n—l 


f{z) 


z=iA 


[n 


1)! dz‘ 


n—l 


ixy 


z=—iA 


= - (p dz- 

27ri 


f{z) 


(13) 

The right hand side of this equation is the width calculated to leading order in 1/M in the 
effective theory, Fq; averaged over masses with a particular weight. We have introduced a 
parameter n, the power of the denominator in ([l3|) , to guarantee vanishing of the integral 
on the circle at inhnity. It needs to be adjusted depending on the number of spacetime 
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dimensions. Defining 


( 14 ) 


/ OO 

dxw{x)f{x) 

-OO 

with the weight function defined by 

^ (n-1)! 1 / 2A^ 

(2u - 3)!! 27rA - M)2 + A2 

and recalling that the width vanishes when ms < rriu + we have obtained 

(f fMt\ = 

' ” (2ii-3)!! + 

This is of the form of our main result, but is not quite complete. It gives a peculiar average 
over heavy quark mass of the decay width to leading order in 1/M in terms of the off- 
shell effective theory Green function with complex momentum. The right hand side can 
be evaluated using an operator product expansion if A is large enough. As in the case of 
semileptonic decays^ the leading operator in the expansion is of the form h^Th^ and has a 
normalized matrix element in the state H{v). Since its coefficient is computed perturbatively, 
one has, to leading order in 1/M, 

(f„(Af)> = (TgiM)) (17) 

where rQ(M) is the perturbative width of the heavy quark. 

The next term in the operator product expansion of T{iA) involves an operator with one 
derivative, which has vanishing expectation value in the state H{v). Therefore, the 

leading correction to Eq. (|T7|) from the OPE is of order 1/M^. However, the Hamiltonian 
itself has an expansion in 1/M so we have yet to establish the validity of quark-hadron duality 
for the decay width to order 1/M^. Moreover, the relation between the weak Hamiltonian 
and its representation in the HQET involves Wilson coefficients that have explicit mass 
dependence. Both types of correction spoil the invariance M —> M -f- 6M, Q ^ Q — 6M. 
We address these issues next. 

The effective Hamiltonian Ti has an HQET expansion in powers of 1/M, 

Ti. = CqHq -\- + ■ ■ ■ (18) 

A sum over several possible operators at each order in 1/M is implicit. The coefficients 
Ck are functions of M/p, where p is a renormalization point which we chose to be a fixed 


Jn—1 


+ 


T(z) 


z=iA 


dz 


n—1 


iAy 


( 16 ) 


z=-iA, 
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FIG. 3: Contour for the integral in the complex Q plane leading to Eq. p^). The thick lines on 
the real axis denote the cuts in the Green function and coefficient function Ck- 


number, large enough that the coefficients can be computed perturbatively. We do not set 
H = M since this would introduce additional M dependence into the operators (which are 
also renormalized at the scale /i). There is a corresponding expansion of the Green function 
in Eq. (j^), 

T = Cofo + + ■ ■ ■ (19) 

and of the width, 

r = fo + fi + ■ ■ ■ ( 20 ) 

Consider the individual averages 


{n(M)> ^ I dQ m{Q)tt(Q) 


( 21 ) 


where the weight function w is given in Eq. (|T^. In order to use a dispersion relation like 
in dl^) we note that the explicit inverse powers of mass give poles at ^ = —M and the 
Wilson coefficients, with typical ln(M) behavior, give cuts extending from z = —M to —oo. 
Using the contour in Fig. |^, which excludes these cut and pole, we are led to consider 


1 / Ck{z + M)Tii.(z) _ 

2Tri f {z M)^{z‘^ + A^)"- 

1 Ck{z + M)fk{z) 

{n — 1)! dz'^~^ {z + MY{z + iA)"- 


+ 


( 22 ) 


1 Ci{z + M)n(z) 

[n — 1)! dz'^~^ (z + MY{z — iA)^ 
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On the other hand, the integral can be written as a sum of two terms, namely, the width 
average we want and the integral over a contour below and above the cut on the negative 
real axis. The latter is suppressed by powers of M. To estimate it we note that since the 
Green function T is analytic in this region we may simply replace it by a power of the mass 
given by dimensional analysis, T{z) ~ (z + My where p = 5 in four dimensions {p = 2D — 3 
in the general case of D dimensions). Also, we may take the Wilson coefficient to be a simple 
log for this estimate. Then the integral around the cut is 


1 

27ri 


-—: / dz 


ln{z + M){z + My-’^ 

{z^TA^ 

1 \ \n{z + M){z + My~^ 


{n — 1)\dz^~^ {z + iAy 


+ (A ^ -A) 


z=iA 


For comparison, a similar estimate of the average (rfc(M)) using rfc(a;) ^ ^ gives 


dx-, 


X‘ 


.p—k 


2>r (2n fl + 0 (h) 


(23) 


(24) 


Jo [{x — My + A'^y (2A)2"-“^ [(n — 1)!]^ 

Thus, we can express the width average in terms of the off-shell Green function and deriva¬ 
tives in Eq. up to corrections suppressed by (A/M)^"’“^ ln(M). 

As in the case of semileptonic decays^, we can now calculate the width average (r(M)) 
by computing the off-shell Green functions using an OPE. But we can do better: we can show 
that there are no corrections of order 1/M. As in the case of semileptonic decays, the OPE 
is an expansion in operators of the form h^VD^h^ where H is a covariant derivative. The 
central observations are that the operator has vanishing expectation value0 and 

the expansion for Tk starts at order h^TD^h^. The latter statement is non-trivial. Gonsider 
the case of semileptonic decays. In the original work of Ref. 0 this question is sidestepped 
by doing a simultaneous expansion in large momentum transfer (the OPE) and the HQET. 
But, following Ref. E3, one could hrst express the Green function in the HQET and only 
then do the OPE. Although Ref. 0 does not consider the effect of subleading operators, 
it is clear that the two approaches yield the same result only if the OPE of products of 
subleading operators starts at the corresponding order in Glearly, a similar indirect 

argument applies here, but we know of no direct proof of the statement. 

Our main result is then 

(25) 


(r(M)) = {rg(M)) (1 + O(j^) 
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It should be noted that the corrections that have been omitted are parametrically small at 
large M, but can be quantitatively large, depending on the values of M, n and A. 

In the’t Hooft model our result is in agreement with the empirical observations of Ref. 
where the averaged widths agree to order 1/M^ while the un-averaged ones agree at best 
to order 1/M. The weight functions used for averages in Ref. [Q were Gaussian, w{x) ~ 
exp(—(x —M)^/A^). We have checked that the results still hold for the weight function in 
Eq. (^). It is interesting that substantial duality violation is found if the power n = 1 is used. 
In the realistic case of QCD in four dimensions one is left with the very realistic possibility 
that the physical hadronic width of a heavy meson exhibits oscillations of magnitude 1/M 
about the partonic width which are erased out when performing unphysical mass averages. 


Some evidence for this was presented in Ref. []^ where it was observed that the 6-quark 
width agrees better with experimental hadronic widths if the quark mass is replaced by 
the B or Af, masses, respectively, and in Ref. [113 which argues that the — Dg lifetime 


difference is also primarily a phase space effect. In a similar vein. Ref. [lAj shows how 1/M 
violations to local, but not global, duality may occur in R-meson correlations. 

To summarize, we have shown that the hadronic width of a heavy meson averaged over 
the heavy quark mass as in Eq. ([T^) is correctly given by the corresponding average of a 
perturbative heavy quark width up to corrections of order 1/M^. This result can be applied 
to the decay widths of heavy mesons in the ’t Hooft model, and explains the numerical 
observations of Ref. |^. The result, however, is not of direct phenomenological signihcance 
since it is impossible to perform mass averages of observed decay widths of B mesons. 
However, our result adds to the body of evidence that heavy meson widths cannot be 
reliably computed using perturbation theory, at least not with a precision of order 1/M^. 
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